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Using a sample of 106 million ^(3686) events collected with the BESIII detector at the BEPCII 
storage ring, we have made the first measurement of the Ml transition between the radially excited 
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charmonium S-wave spin-triplet and the radially excited 5-wave spin-singlet states: i/;(3686) — > 
jr] c (2S). Analyses of the processes V(3686) -> 7??c(2S') with r/ c (2S) -> A"giC ± 7r =F and K + K~tv° 
gave an r/ c (2S) signal with a statistical significance of greater than 10 standard deviations under a 
wide range of assumptions about the signal and background properties. The data are used to ob- 
tain measurements of the rj c {2S) mass (M(rj c (2S)) = 3637.6 ± 2.9 sta t ± l-6 sys MeV/c 2 ), width 
(T(nc(2S)) = 16.9 ± 6.4s ta t ± 4.8 sys MeV), and the product branching fraction (23(0(3686) -> 
7?7 C (2S)) x B(r] c (2S) -> KKn) = (1.30 ± 0.20 sta t ± 0.30 sys ) x 10~ 5 ). Combining our result with a 
BaBar measurement of B(rj c (2S) — > KKn), we find the branching fraction of the Ml transition to 
be 23(0(3686) -> -JVc(2S)) = (6.8 ± l.l stat ± 4.5 sys ) x 10" 4 . 

PACS numbers: 13.25.Gv, 13.20.Gd, 14.40.Pq 



Recent discoveries of charmonium and charmonium- 
like states above the open-charm production threshold 
have generated great interest. Intensive efforts to in- 
corporate these states into the quark-model picture of 
hadrons have led to the development of models to ex- 
plain some or all of the new states The charmo- 
nium states below the open-charm production threshold 
are relatively well understood, with the notable excep- 
tion of the spin singlets. These include the P-wave state 
h c and the iS-wave ground state r\ c and its first radial 
excitation rj c (2S) 0]. These are experimentally challeng- 
ing because of the low production rates and spin-parity 
quantum numbers that are inaccessible in direct e + e~ 
annihilations. 

The r] c (2S) was first observed by the Belle collab- 
oration in the process B ± —> K ± r) c (2S), ij c (2S) —> 
KgK^irT It was confirmed in the two-photon 

production of KgK^n^ [1, [Bj], and in the double- 
charmonium production process e + e~ — > J/ipcc @, @]- 
Combining the world-average values Q with the most re- 
cent results from Belle and BaBar on two-photon fusion 
into hadronic final states other than KgK^it^ one 
obtains updated averages of the rj c (2S) mass and width 
of 3637.7±1.3 MeV/c 2 and 10.4±4.2 MeV, respectively. 

The production of the rj c (2S) through a radiative tran- 
sition from the "0(3686) requires a charmed-quark spin- 
flip and, thus, proceeds via a magnetic dipole (Ml) 
transition. The branching fraction has been calcu- 
lated by many authors, with predictions in the range 
23(0(3686) -> 77? c (25)) = (0.1 - 6.2) x 1CT 4 0. A 
recent calculation [llj that includes contributions from 
loops containing meson pairs finds a strong cancellation 
that results in a partial width of (0.08 ± 0.03) keV and 
a branching fraction of (2.6 ± 1.0) x 10 -4 ; while a calcu- 
lation using the light-front quark model and a 2S state 
harmonic oscillator wave function to present the 2S char- 
monium state gives a transition rate of 3.9 x 10~ 4 fl2j |. 
Experimentally, this transition has been searched for by 
Crystal Ball [lj], BES 0], CLEO [lland most recently 
by BESIII through r] c (2S) -> [lf|. No convincing 
signal was observed in any of these searches. 

In this Letter, we report the first observation of 
0(3686) -> mc(2S), with r] c (2S) -> KgK^ and 
K + K~ir°. The data sample for this analysis consists of 



an integrated luminosity of 156 pb _1 (106 million events) 
produced at the peak of the 0(3686) resonance [ItJ and 
collected in the BESIII detector [18j. An additional 
42 pb _1 of data were collected at a center-of-mass energy 
of ^—3.65 GeV to determine non-resonant continuum 
background contributions. 

The BESIII detector, described in detail in Ref. fl8|], 
has an effective geometrical acceptance of 93% of 47t. 
A small-cell, helium-based main drift chamber (MDC) 
in a 1-T magnetic field provides a charged-particle mo- 
mentum resolution of 0.5% at 1 GeV/c, and specific- 
ionization (dE/dx) measurements for particle identifica- 
tion with a resolution better than 6% for electrons from 
Bhabha scattering. The cesium iodide electromagnetic 
calorimeter (EMC) measures photon energies with reso- 
lutions at 1.0 GeV of 2.5% and 5% in the detector's barrel 
(|cos#| < 0.8, where 9 is the polar angle with respect to 
the e + direction) and endcaps (0.86 < |cos6>| < 0.92) re- 
gions, respectively. Additional particle identification is 
provided by a timc-of-flight system (TOF) with a time 
resolution of 80 ps (110 ps) for the barrel (endcaps). 

Reconstructed charged tracks other than daughters of 
Kg candidates are required to pass within 1 cm of the 
e + e~ annihilation interaction point (IP) transverse to the 
beam line and within 10 cm of the IP along the beam 
axis. Each track is required to have a good-quality fit 
and to satisfy the condition |cos#| < 0.93. Charged- 
particle identification (PID) is based on combining the 
dE/ dx and TOF information to construct a Xpid 

(i). The 

values XpidW an< ^ the corresponding confidence levels 
ProbpiD(i) are calculated for each charged track for each 
particle hypothesis i (pion, kaon or proton). 

A neutral cluster in the EMC must satisfy fiducial and 
shower-quality requirements to be accepted as a good 
photon candidate. Showers must have a minimum energy 
of 25 MeV and be detected in either the barrel or endcap 
regions, as previously defined. EMC timing requirements 
are used to suppress noise and energy deposits unrelated 
to the event. 

In selecting ^KgK^n^ ( / -fK + K~Tr°) events, the decay 
signal K s — > tt + tt~ (tt° — > 77) is used to tag the K s 
(tt°). Candidate events must therefore have exactly four 
(two) charged tracks with zero net charge and at least one 
(three) good photon(s) for the 7^g^ ± 7r =F (-/K+K'-tt ) 
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decay mode. 

Kg candidates are selected with secondary- vertex fits 
to all pairs of oppositely charged tracks in the event, 
assuming pion masses. The combination with the best 
fit quality is chosen and the event is kept for further 
analysis if the invariant mass is within 7 MeV/c 2 of the 
expected Kg mass, and the secondary vertex is at least 
0.5 cm from the interaction point. To suppress ^jK s K s 
events, the remaining tracks are required not to form a 
good K s candidate. The fitted K s information is used 
as input for the subsequent kinematic fit of the complete 
event. 

The r yK s K ± Tr T candidates are then subjected to a 
four-constraint (4C) kinematic fit, with the constraints 
provided by four-momentum conservation. The dis- 
crimination of charge-conjugate channels (KgK + ir~ or 
KgK~ir + ) and the selection of the best photon among 
multiple candidates arc achieved by minimizing x 2 = 

xlc +Xpid( k ) + Xpid( 7F )> where xlc is the chi-square of 
the 4C kinematic fit. Events with x\c < 50 are accepted 
as 7 AgA ± 7r =F candidates. For jK + K~n° candidates, 
both charged tracks must satisfy the criterion that the 
kaon- hypothesis probability ProbpiD(A) is larger than 
both 0.001 and the probability of any other hypothesis. 
A five-constraint (5C) kinematic fit, with the ir° mass 
as the additional constraint, is used to select the best 
transition photon and the ir a —> 77 combination. Events 
with Xsc < 30 are accepted as jK + K~ir° candidates. 

We use the program lundcrm [3] to generate inclu- 
sive Monte Carlo (MC) events for background studies. 
The signal is generated with the expected angular dis- 
tribution for -0(3686) — > jr) c (2S), and the subsequent 
i] c (2S) — > KgK ± n T and K + K~~ir° decays are generated 
according to phase space. The detector response is sim- 
ulated with a GEANT4-based package [2(| that has been 
tuned to match the performance of the detector compo- 
nents. 

The 0(3686) -> 7?7 C (2S*), 77 C (2S*) -> K^K^ 
(K + K~tt°) signal suffers significantly from background 
contributions from leptonic J/ip decays and J/ij) — > 
K+K- in 0(3686) -> 7T+7T-J/0, and 0(3686) -> 77 J/0 
with rj — > 7T + 7r~7r° (77). For the 7 AijA ± 7r T channel, 
these background contributions are suppressed by re- 
quiring that the recoil mass of all ir + ir~ pairs be less 
than 3.05 GeV/c 2 . For the jK+K~tt° channel, this type 
of contamination is removed by requiring that the in- 
variant mass of the two charged tracks, assuming they 
arc muons, to be less than 2.9 GcV/c 2 . The remain- 
ing dominant background sources are (1) 0(3686) — !> 
KgK ± ir T (K + K~ir°) events with a fake photon can- 
didate; (2) events with the same final states including 
-K'gA' ± 7r =F 7 ISR / FSR (if + A"-7r°7 ISR/FSR ) with the pho- 
ton from initial- or final-state radiation (ISR, FSR) and 
0(3686) — > ujK + K~ with us —> -fir ; and (3) events with 
an extra photon, primarily from 0(3686) — > n° K^K^ir^ 



(tt°K + K~tt ) with 7T° — > 77. MC studies demonstrate 
that contributions from all other known processes are 
negligible. 

The events in the first category, with a fake photon 
incorporated into the kinematic fit, produce a peak in 
the A^A ± 7r =F (K + K~tt°) mass spectrum close to the 
expected r] c (2S) mass, with a sharp cutoff due to the 
25-MeV photon-energy threshold. 

Because the fake photon adds no information to the fit, 
its inclusion distorts the mass measurement. We there- 
fore determine the mass from a modified kinematic fit 
in which the the magnitude of the photon momentum 
is allowed to freely float (3C for r yKgK ± ir T and 4C for 
r yK + K~ / K°). In the case of a fake photon, the momentum 
tends to zero, which improves the background separation 
with minimal distortion of the signal line shape fl(| . 

Background contributions from 0(3686) — > KgK ± TT T 
{K+K-ir°) and 0(3686) -> AgA ± 7 rT 7FSR 
(A + A _ 7r 7F s R ) are estimated with MC distribu- 
tions for those processes normalized according to a 
previous measurement of the branching ratios [2l| . 
FSR is simulated in our MC with PHOTOS and 
the FSR contribution is scaled by the ratio of the 
FSR fractions in data and MC for a control sample 
of 0(3686) — > 7Xcj (J = or 1) events. For this 
study the \cj is selected in three final states with or 
without an extra FSR photon, namely A^A ± 7r^(7 F g R ), 
7r + 7r~7r + 7r~ (7 F sb), and 7r + 7r~A' + A'~(7 F s R ), as de- 
scribed in Ref. [l6| . Background contributions from 
the continuum process e + e~ — > 7 * — > A^A ± 7r =F (7 F s R ) 
(A + A~7r ( 7F s R )) and the ISR process e + e~ — s- 
7*7isr — > A'gA ± 7r T 7 is R (A + A~7T 7 is R ) are estimated 
with data collected at y/s = 3.65 GcV corrected for 
differences in the integrated luminosity and the cross 
section, and with particle momenta and energies scaled 
to account for the beam-energy difference. MC simula- 
tions show that the AgAT ± 7r =F (K + K~tt°) mass spectra 
are similar for FSR and ISR events. Events without 
radiation have the same mass distribution independently 
of originating from a resonant 0(3686) decay or from 
the non-resonant continuum production. Thus, the 
background shapes from KgK^w^ (A + A~7r°) and 
A's A ' ±7rT 7iSR./FSR (A+A"7r°7 ISR/FSR ) are described by 
the sum of the MC-simulated K^K ± n T (A+A^tt ) and 
AgA ± 7r T 7 F s R (A + A~7r°7 F s R ) invariant-mass shapes, 
with the proportions fixed according to the procedure 
described above. The shapes of backgroung mass 
distributions from 0(3686) — » ojK + K~ with oj — > 7 7r° 
are parameterized with a double-Gaussian function, and 
its level is measured with the same data sample and 
fixed in the final fit. 

The third type of background, that with an extra pho- 
ton, ir°K^K ± ir T (7r°A + A'~7r°), is measured with data 
and normalized according to the simulated contamina- 
tion rate. It contributes a smooth component around 
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the Xcj (J = 1,2) mass region with a small tail in the 
ij c (2S) signal region that is described by a Novosibirsk 
function |23[ (Gaussian function) for the tt°K s K ± tt t 
(ir°K + K~~ir°) background. The shape and size of this 
background is fixed in the fit. 

The mass spectra for the KgK ± TT T and K + K~n° 
channels are fitted simultaneously to extract the yield, 
mass and width of r] c (2S). To better determine the back- 
ground and mass resolution from the data, the mass spec- 
tra are fitted over a range (3.46-3.71 GeV/c 2 ) that in- 
cludes the Xci arid Xc2 resonances as well as the rj c (2S) 
signal. The final mass spectra and the likelihood fit re- 
sults are shown in Fig. [T] Each fitting function includes 
four components, namely the r] c (2S), Xci, Xc2, and the 
summed background described above. Line shapes for 
the Xci an d Xc2 are obtained from MC simulations and 
convolved with Gaussian functions to accommodate for 
the mass-scale and resolution differences from data. For 
both modes, the Xci and Xc2 widths are fixed to the PDG 
values 0. Based on MC studies, the mass shift and res- 
olution for the resonances are found to vary linearly as 
a function of the KgK ± n T (K + K~tt°) invariant mass. 
These parameters are extrapolated from the Xci and Xc2 
to the r] c (2S). 

The line shape for the r] c (2S) produced in the Ml 
transition of the 0(3686) is assumed to have the form 
x BW(m) x fd(E-y) x e(m)) <x> G(Sm, a), where m 
is the invariant mass of the KgK+n^ or K + K~tt°, 
e i = ( m v-(3686) ~ to2 )/ 2to v>(3686) is the energy of the 
transition photon in the rest frame of "0(3686), BW(m) 
is the Breit-Wigner function for r] c (2S), fd(E 7 ) is a func- 
tion that damps the diverging tail originating from the 
E^ dependence, e(m) is the mass-dependent efficiency 
function determined by a full simulation of the signal, 
and G(Sm, a) is a Gaussian function describing the mass 
shift and the detector resolution. For the damping func- 
tion we use a functional form introduced by the KEDR 
collaboration H|: fd = E^/[E-yE Q + (E 1 - E a )] 2 , where 
E n = {™% ( 368 6 ) - m ^ ( 2S))/ 2m V'(3686) is the peak energy 
of the transition photon. To assess the sensitivity of our 
results to the choice of this function, we also consider an 

alternative form used by CLEO [25|: fd = cxp(— ^j-), 
with CLEO's fitted value of /3 = (65.0 ± 2.5) MeV. 

The fit shown in Fig. Q] has a x 2 of 72 for 79 degrees 
of freedom. The results for the yields of r] c (2S) events 
are 81 ± 14 for the K^K^-k^ channel and 46 ± 11 for 
the K + K~tt° channel. Consistent yields are found for 
separate fits to the two channels [26j. The KgK^ir^ 
channel determines primarily the precision for the rj c (2S) 
mass and width measurements in the simultaneous fit 
with the results M Va{2 s) = 3637.6 ± 2.9 MeV/c 2 and 
rr; c (2S) = 16-9 ± 6.4 MeV, respectively. The combined 
statistical significance of the signal in the two modes is 
ll.lcr, which is obtained by comparing the likelihoods of 
the fits with and without the r/ c (2S) signal. The robust- 



ness of this result was tested by considering variations 
of the resonant line shapes, background assumptions and 
other systematic effects. In all the cases, the statistical 
significance is found to be larger than 10.2o\ 

Combining the observed number of signal events with 
the efficiencies of 25.6% and 20.2% for the K° s K ± tt+ and 
K + K~tt° final states, respectively, from full simulations 
of the signal with the measured rj c (2S) mass and width, 
we find the product branching fractions £>(?/;(3686) — >• 
7?7 C (2S')) x B{r) c {2S) -> K^K^) = (4.31±0.75) x 10~ 6 , 
and £(0(3686) -> -/r] c {2S)) x B{r/ c (2S) -> K+K-jt ) = 
(2.17 ± 0.52) x 10~ 6 , where the errors are statistical 
only. The ratio of the branching fractions agrees well 
with the isospin symmetry expectation of 2:1 between 
KgK ± n T and K + K~ir°. The product branching frac- 
tion for 0(3686) -> jr) c (2S), r] c (2S) -> KKir can be 
obtained by doubling the sum of the KgK ± ir T and 
K + K~~ir° branching fractions to obtain i3(^(3686) — > 
-yrj c {2S)) x B{rj c {2S) -> KKir) = (1.30 ±0.20) x 10~ 5 , 
where the error takes into account the correlation be- 
tween the two measured branching fractions from the si- 
multaneous fit. 

The systematic uncertainties in the branching frac- 
tion, ?7 C (25") mass and rj c {2S) width measurements are 
summarized in Table [T] The uncertainties due to the 
choice of the background shape, the damping func- 
tion, the fitting range and the linear extrapolated mass 
shift for ?7 C (2S') are common among the three mea- 
surements and arc determined together. The sys- 
tematic errors in the mass and width due to the 
- ftr s A ' ±7rT (7iSR/FSR.) ( a+a_ ^°(7isr/fsr)) background 
shape are evaluated by changing the relative ratio of 
the K'gK ± -K^ (K + K~tt°) background events with and 
without radiation. The uncertainties from ir° KgK ± 7T T 
(tt°K + K~tt ) background shape are estimated by chang- 
ing the function parameterizing the measured mass spec- 
trum. The uncertainty due to the choice of damping func- 
tion is estimated from the difference between results ob- 
tained with the default (KEDR) and alternative (CLEO) 
functional forms. The uncertainties due to the choice 
of fitting range are estimated by taking the largest dif- 
ferences between results found with the standard fitting 
range and those obtained using alternative ranges. The 
uncertainties from the linear extrapolation of the mass 
shifts from Xci and Xc2 to r] c (2S) are estimated from the 
maximum changes in the fitting results obtained by vary- 
ing the mass shifts within their errors. 

The branching fraction measurement is affected by 
additional effects that enter through the yield determi- 
nation, including those associated with charged-particle 
tracking, photon reconstruction, particle identification, 
Kg reconstruction, and kinematic fitting (% 2 require- 
ment) , all of which are estimated with control samples in 
the data j2?| . The effect of the uncertainty in the dynam- 
ics of the decay rj c (2S) -> K s K ± tt^ (K+R-tt ), which 
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FIG. 1: The invariant-mass spectrum for K%K ± n^ 1 (left panel), K + K ir° (right panel), and the simultaneous likelihood fit to 
the three resonances and combined background sources as described in the text. 



TABLE I: The absolute systematic uncertainties in the r] c (2S) 
mass (in MeV/c 2 ), width (in MeV) and the relative sys- 
tematic error (in %) in BB, the product branching fraction 
23(^(3686) -> 7^(25)) x B(rj c (2S) -> KKir), measurements. 



Source 


Mass 


Width 


BB 


Background shape 


1.3 


2.6 


9.9 


Damping function 


0.7 


4.0 


19.6 


Fitting range 


0.1 


0.4 


1.3 


Mass shift 


0.6 


0.2 


0.4 


Tracking 






4.0 


Photon reconstruction 






1.3 


Particle identification 






1.3 


Kg reconstruction 






2.3 


Kinematic fitting 






3.9 


r) c {2S) decay dynamics 






1.5 


Number of 0(3686) events 






4.0 


Total 


1.6 


4.8 


23.3 



is treated as phase space in our default signal MC, is esti- 
mated with an alternative MC replicating the Dalitz dis- 
tribution of rj c (2S) — > KgK^-K 11 decay recently measured 
by the Belle collaboration A 0.8% (3.0%) relative 

difference in the efficiency was found between the default 
and alternative MC samples for K^K^tt^ (K + K~tt q ), 
leading to a 1.5% difference in the total branching ratio, 
which we take as a systematic error. Finally, there is 
an overall 4% uncertainty in the branching fraction as- 
sociated with the determination of the total number of 
-0(3686) events in our data sample fl7| . 

We assume that all the sources of systematic uncertain- 
ties are independent and combine them in quadrature to 
obtain the overall systematic uncertainties given in Ta- 
ble HI The total systematic uncertainties on the mass and 
width measurements are 1.6 MeV/c 2 and 4.8 MeV, re- 
spectively; the total relative systematic uncertainty on 
the product branching fraction £>(-0(3686) — > 7^(25')) x 



B(i] c (2S) —> KKir) is 23.3%. Using the measurement of 
B(r? c (2S) -> KKir) = (1.9 ± 0.4 ± 1.1)% from the BaBar 
experiment [jyl ]. we find an Ml-transition branching frac- 
tion of £(^(3686) -> ir] c (2S)) = (6.8 ±1.1 ±4.5) x lO" 4 , 
where the systematic error is dominated by that of the 
BaBar result. 

In summary, we report the first observation of the 
Ml transition 0(3686) — > 7?y c (25') through the decay 
processes V(3686) -> yK^K*^ and jK+R-tt . We 
measure the mass of the r/ c (2S) to be 3637.6 ± 2.9 ± 
1.6 MeV/c 2 , the width 16.9 ± 6.4 ± 4.8 MeV, and the 
product branching fractions 23(0(3686) — > -fr] c (2S)) x 
B(r] c (2S) -> KKir) = (1.30 ± 0.20 ± 0.30) x 10" 5 , where 
the quoted uncertainties are statistical and systematic, 
respectively The main systematic limitations to these 
measurements arise from the choice of the functional 
form for the damping factor in the rj c (2S) line shape 
and from uncertainty in the choice of the background 
line shapes. Our results are consistent with previously 
published values and limits, and the branching-fraction 
measurement of the Ml transition 0(3686) — > jr] c (2S) 
of (6.8 ± 1.1 ± 4.5) x 10~ 4 agrees with theoretical cal- 
culations and naive estimates based on the J/0 — > "frj c 
transition (l5| . 
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